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Abstract — Ship autopilots are usually designed based on PID 
controller because of the simplicity and the ease of construct. 
However its performances in various environmental conditions 
are not as good as desired. To improve the ship performances, a 
linear state space feedback autopilot is designed. The paper 
presents the utility of this autopilot to stabilize the ship, shaping 
its yaw response and force the system to follow a desired path. 
The state space model of the 4DOF ship motion is presented and 
the method of the autopilot designing is completely described. 
Furthermore numerical simulations using real ship parameters 
are provided to illustrate the effectiveness of the proposed 
methodology. For the ship path following, a comparison between 
PID and feedback gain controllers is given to show the robustness 
of each controller in presence of wave disturbances. 
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I. Introduction 

Marine vehicles are designed to operate with adequate 
reliability and economy. And in order to accomplish this, it is 
essential to control the motion of the ship. Despite the external 
disturbances (waves, wind, currents), the control task consists 
in making the ship to follow, as closely as possible, a desired 
trajectory, which can be defined in terms of the ship's position, 
velocity and acceleration [1]. 

Over the last two decades, the ship control had received 
considerable research efforts and had undergone a period of 
significant progress, especially for roll and yaw motions. This 
was the work of different groups of research which targeted the 
development in different countries as: The Netherlands, 
Denmark, Sweden, the USA, Japan and China. 

Conventional ships which are moving with constant 
forward speed, the steering machine is the only system to 
control the ship's motion. The principal function of the steering 
machine is to move the rudder to a desired angle when 
demanded by the autopilot system or by the helmsman. The 
main purpose (or the helmsman) of the autopilot is to control 
the heading of the ship affected by disturbances [2]. 

In actual ship building field, a considerable number of the 
vessels traveling oceans of the world are only equipped with 
rudders [3]. In this paper, it is assumed that the ship is moving 
at a constant forward speed, and only the rudder is controllable. 
The ship control problem is how to make the ship keep the 
desired trajectory by controlling the rudder's actuator. 

To design an autopilot for the ship's motion is always a 
challenging problem. The motion is influenced by 
unpredictable environmental disturbances such as waves, winds, 
currents, change of depth under keel, etc. as well as ship sailing 
conditions such as speed, loading condition, trim, etc. To 
design an autopilot capable of taking into account all these 
factors is nowadays possible by different control design 



methods [4], The autopilot takes into account the yaw-angle 
variations and it generates rudder commands to reject the 
disturbances according to the heading control law. 

Many control methods are used to improve the ship 
performances especially against the disturbances, but the 
implementation of these methods are still difficult for 
conventional systems [5, 6], However, state-space feedback 
gain, as one of the simplest methods, requires only a linear 
state-space model of the plant to be applied and implemented 
for a controllable system. 

Because the necessity of the system model, existing 
mathematical ship models were developed to incorporate the 
use of all kind of controllers that could be commanded in a 
strategic manner. The models are still complex to use in control; 
many studies tried some modifications and simplifications to 
make it useful and easy to be controlled. Adopting these 
models, some theoretical study results have been achieved on 
the ship control motion [6]. 

In this work, a linear model for underactuated ship motion 
is presented in horizontal plane (4 Degree Of Freedom). This 
model is developed to state space representation that allows 
applying the feedback controller to shape ship motion as 
desired. To design the feedback controller, we will focus on the 
pole placement method. The procedure consists of two main 
steps: (i) proving the controllability of the system then writing 
its canonical form; (ii) finding the state-feedback gain by fixing 
the polynomial of the desired response. The principal role of 
the controller is to improve the ship performances in terms of 
effectiveness and robustness whatever the sailing conditions. 
The simplicity of the controller, comparatively with the modem 
and intelligent theory of control, can be deduced during the step 
of the feedback controller design. 

The rest of the paper is organized as follows. System of 
interest and problem formulation is presented in the next 
section. Section 3 devotes to state-space modeling and section 
4 describes the feedback control. Design procedure is described 
in section 5 while results and discussion are given in Section 6. 
Section 7 concludes the paper. 

II. System of Interest and Problem Formulation 
It is known that to fully represent a rigid body motion in 

space requires a Six Degree Of Freedom (6DOF) approach. 

This is the same case for the ship as a rigid body moving in sea. 

Fig. 1 represents the Standard Notation followed by the 
Society of Naval Architects and Marine Engineers notation 
convention (SNAME. 1950) applied on the ship motion 
(6DOF). 

This notation employs a combination of two separate right- 
handed orthogonal coordinate systems: one attached to the 
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earth and a second attached to the body of the ship. The 
TABLE I describe the different parameters of the Fig. 1. 

The equations of motion describing the dynamics of a ship 
are readily obtained fromNewton's law [5]: 



d(mV) 



E Forces 
dt 

Moments - 



(1) 



dt 



Where: m : mass of the body. V : translational velocity. I : 
Inertial moment. a> : rotational velocity. 



Roll 

p,K 




Heave 
V z, w, z 



q t m Sway 
y, v, Y 



Figure 1 Standardnotation and sign conventions for ship motion description 

(SNAME, 1950) 

TABLE I Coordinate Frame Vessel parameters 



Paramete rs 


Description 


x,y,z 


Distance along body-fixed system axesfm] 


u ,v,w 


Translational velocity in body-fixed frame [ml s] 


Xf,Yf,Zf 


Force components along body axes [kg] 


p,q,r 


Rotational velocity in body-fixed frame [rad I s] 


K,M,N 


Moment components along body axes [kg .m] 


¥ 


Yaw angle [rad] 


e 


Pitch angle [rad] 


<t> 


Roll angle [rad] 


8 


Rudder angle [rad] 



Applying (1) on the ship, the motion in earth-fixed frame is 
described as followed: 

f x (u, u, v, v, w, w, p, p, q, q, r, r, 8) = Xf 

f y (u, u, v, v, w, w, p, p, q, q, r, r, 8) = Yf 

f z (u, u, v, v, w, w, p, p, q, q, r, r, 8) = Zf 

f^(u,ii,v, v, w, w, p, p, q,q,r,r,8) = K (2) 

f g (u, u, v, v, w, w, p, p, q, q, r, r,8) = M 

f v (u, u, v, v, w, w, p, p, q, q, r, r,8) = N 

However a 3DOF plane motion is usually considered 
adequate for ship maneuvering study (surge, sway and yaw) [7], 
and for high speed vessels, roll mode is not negligible. Hence a 
4DOF description that includes surge, sway, yaw and roll 
modes is needed [8], (Fig. 2). 

Upon linearization with respect to straight line motion with 
forward speed u , the surge equation is decoupled and the 

following linear coupled sway-yaw-roll equations follow 
immediately [8]. 



For simplicity we assume that the only external force and 
moment are caused by a single rudder angle noted bye?, where: 

[Yf N K Of =[Y S N s X, Of (3) 
Referring to Fig. 2, we have: 



m(v + u r) = Y i v + Y v v + Yj + Y r r + Y p p + Y p p + Y^ + Y lj S 

if/ = r 

IJ = K v v + K x y + K r r + Kj + K p p + K p p-mgGMif> + K s S 

<f) = p 



(4) 



o ( 



••-' 8 




Figure 2 surge-sway -yaw -roll motion coordinate system 

Where: 

Y v ,N v ,K. : indicate the hydrodynamic coefficients. 

For instance, Y indicates the derivative of sway force Y to the 
sway acceleration v evaluated at the reference condition. 

m : is the mass of the ship [faj.s 2 Irn] ■ 

I I : are the moments of inertia about x-axis and z-axis 
respectively [kg s 2 .m]. 

g Acceleration of gravity [ml s 2 ]. 

GM : is the metacentric height [m], which indicates the 
restoring capability of a ship in rolling motion [9]. 
From (4) we rewrite: 

v = Y'v + Y'v + Yf+Y'r+Y'p + Y'p + Y'd + Y'8 

v v r r pi p r $T d 

f = N'v + N'v+N'f + N r + N'.p + N' p + N'J+N'8 

v v r r pr pi ft d) 

W = r (5) 

p = K' v v + K'v + K'j + Kj + K' p p + K' p p + K'^ + K' s 8 

\j> = p 

The equations in (5) present the linear model of the coupled 
steering and roll dynamics. 
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III. State space Modeling 

State space representation is a mathematical model of a 
physical system as a set of input, output and state variables 
related by first-order differential equations. The state space 
representation provides a convenient and compact way to 
model and analyze systems with multiple inputs and outputs. 

For notational convenience, we define the state vector as: 

x = [v r iff p </>J (6) 

Equation (5) can be written in matrix form as followed: 



v \Y' V Y' Y' 

r N' N' N' N[ 

\jr = 10 

p K' K' K' K' 

r v r p ft 

ij>\ [o 1 
Hence: 
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(8) 



If M 'exists, we have: 



a n a 12 O 

a 21 a 22 O 

OlO 

Q 41 Q 42 ° 

OOO 



Q 14 Q V 

a 24 a 2! 

O 

Q 44 Q 4. 

1 



(9) 



To simplify the system for further analysis, we can 
reorganize the state vector again such that state variables 
associated with the yaw and roll dynamics are separated. 
Moreover, (9) can be rewritten as: 



X 

V 




X, 

f 1 _ 








Separating the dynamics, we can write: 



x =A x 
= A. x ■ 

0V Y 



- A^ x t + b 5 Sway - Yaw dynamics 
A . . x. +b,5 Roll Dynamics 



Neglecting the coupling matrices [2] (A ^ = 0, A. 
implies that: 



(10) 

(11) 
= 0) 

(12) 



Equation (12) is recognized as the Nomoto model. This 
model is widely employed in the ship steering autopilot design. 
Its simplicity and reasonable accuracy in describing small 
rudder angle yaw dynamics make them attractive [9], 

If we consider the yaw angle as the system output (^ = r), 
the full system will be presented in state space model as: 
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(13) 



(14) 



Linearization of the ship nonlinear model about u (service 
speed) implies that we can write the linearized model in 
standard state space form: 

jx(t) = Ax(r) + Bu(r) 

\y(t) = Cx(t) 

One of the important results that we can extract from the 
state space model is the stability of the system. The system is 
stable if and only if the eigenvalues of the matrix A have a 
negative real part. 

We can extract other results as: the controllability, 
observability and identifiability. All of these characteristics can 
be demonstrated for the system. 

IV. Feedback Control 

For the controllable system, the design of state feedback 
control yields desirable closed-loop performance in terms of 
both transient and steady-state response characteristics. If the 
system is not, the state feedback can be used to stabilize the 
closed-loop state equation [10]. 

The state feedback control law features a constant state 
feedback gain matrix K and a new external reference 

input r(t) (Fig. 3). 

Equation (14) represents the open-loop system to be 
controlled. So, the focus is on the application of state feedback 
control laws of the form: 

u(t) = -Kx(t) + r(t) (15) 

With the goal of achieving desired performance 
characteristics for the closed-loop state equation: 



(16) 



jx(t) = (A-BK)x(t) + Br(t) 
\y(t) = Cx(t) 

A fourth order linear state space model representing the 
sway-yaw-roll modes of motion is given in Ref. [11]; we can 
add the equation characterizing the yaw angle {\j/ = r) and the 
full model will be: 



A = 



-0.04 -1.933 -0.0756 
-0.00011 -0.0813 0.001134 

10 
-0.00594 -0.07 -0.059 

1 1 
OlOOO] 



0.1559 
-0.0033 


0.00821 





r(t) + u(t) + 
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x( 


t) 


r *>, 


-j 


t • » ;i 
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lU 







Figure 3 Basic structure of a state-feedback control system 
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V. Design Procedure 

A. Step 1: Controllability 

The controllability is an important property of a control 
system, and it plays a crucial role in many control problems, 
such as stabilization of unstable systems by feedback, or 
optimal control. 

Roughly speaking, the concept of controllability denotes the 
ability to move a system around in its entire configuration 
space using only certain admissible manipulations [12]. 

The first step of designing is to prove that the pair (A, B) is 
controllable. If it is, the vector K exists. 



The controllability matrix is defined as: 

p = [b ab a 2 b a 3 b a 4 b] 



(17) 



The system is controllable if the matrix P has a full-row 
rank (it means: rank(P) = order(A) = 5 ). 

We have: 

0.1559 -0.0129 0.0003 0.0017" 
-0.0033 0.0003 0.0002 

-0.0033 0.0003 0.0002 

0.0082 -0.0015 -0.0096 0.0008 0.0013 

0.1641 -0.0015 -0.0225 0.0011 

rank(P) = 5 , hence the system is controllable and the 
vector K exists. 

B. Step 2: Canonical form 

The second step is to find the matrix T that transforms the 
system in its canonical form. 

The system will be represented as: 



(18) 



J z(t) = (A - BK c )z(t) + B c r(t) 

[y(t) = C e z(t) 
Where: z(t) = Tx(t) 

This gives: 

A =TAT 1 ,B =TB,C =CT 1 and K =KT~ 1 



Note that all the followed results are computed by Matlab 
software. 

For the considered system, it is found that the matrix of 
transformation to the canonical form is: 



-5.08 


-553.41 


-54.82 


-4.12 


-0.29 


-30.83 


-14247.78 


-1578.93 


127.67 


-2.41 


-3159.15 


-145938.58 


-11227.33 


1329.31 


12.48 


-3159.15 


-146365.84 


-14697.75 


1329.31 


-56.37 


-22155.03 


-1021468.36 


-76830.91 


10133.67 


130.18 


Vnd we have: 











A 




10 
10 
10 
1 








1 


0.0007 







0.0206 


B c = 





0.1461 







0.1913 








C c =[-0.0033 0.0002 0.0002 o] 
In next step, we adopt the canonical form of the system and 
our goal will be to find K instead of K . 

C. Step 3: feedback vector 

In this step, the feedback vector K c is expressed in terms of 
the open-loop and closed-loop desired dynamic parameters. 

The open-loop dynamics is characterized by: 

\sl - A | = s 5 + 0.19s 4 + 0.14s 3 + 0.026s 2 + 0.0007s (1 9) 
We denote the feedback gain K as: 

K c = [S S, S 2 S 3 S t ] (20) 

As result, the closed-loop dynamics is characterized by: 



\sl - (A - BK c J = s 5 + (S 4 + 0.19)s 4 + (S 3 + 0.14)s 3 
+ (S 2 + 0.026)s 2 + (S 1 + 0.007)s + S Q 



(21) 



The desired dynamics is expressed by a five-order 
polynomial as followed: 



a(s) 



s 5 +a 4 s 4 + a 3 s 3 +a 2 s 2 +a 1 s 1 + a (22) 



Equalizing the desired dynamics polynomial (22) with 
closed-loop dynamics polynomial (21), we can write the vector 
K as: 



K c = [a (a, - 0.0007) (a 2 - 0.0206) 

(a, -0.1461) (a 4 - 0.1913)] 
Now, it only remains to find the coefficients a. ■ 



(23) 



D. Step 4: Dominant poles approach 

The mathematical formulation of the parameters as: rise 
time t , settling time t , percent overshoot PO etc., is not 
known for high order system ( order > 2 ). Hence, many studies 
tried to approximate the high order system to first or second 
order system if possible. 

The five order polynomial in (22) can be rewritten as: 

a(s) = (s-PJ(s-I[){s-P 2 ){s-P 3 )(s-P i ) (24) 

In this work, the dynamic of yaw ship motion is shaped as 
second order system dynamics (it can be shaped also as first 
order system response). Second order system response is 
usually characterized by two parameters noted E and a> n ■ 

We consider only P o and P of (24) and we neglect P 2 , 
P 3 and P 4 . 

We write: 
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■P 1= -^„+0„jVl-^ 2 (25) 

P 2 =P 3 =P 4 =10*Re(P ) 

Note that: 

j = v-1 . Re : real part of P o 

t, : dimensionless damping ratio. 

ojjrad.s 1 ] : fixed undamped natural frequency. 

For dominant poles approach, we take P 2 , P 3 and P 4 ten 
times P (or more) to neglect their dynamics and the desired 
behavior will obey only to two poles: P o andP t . 

E. Step 5: Shpaing the dyanmics 

The last step is to specify desired closed-loop system 
behavior via (P ,P 1 ) selection; it is called: shaping the 

dynamic response. This choice fixes directly the values of 
feedback gain vector K . 

It consists to fix the couple (g t a> B ) which characterizes the 
desired response. By fixing values of £ and a> , it leads to fix 
the closed-loop characteristics of the transient response such as: 
rise time t R , peak time t p , percent overshoot PO , and settling 
time t s . 

For the ship yaw motion, we have taken the followed values: 

\Z = °- 8 



ft> 



0.02 



(26) 



Where: 



2.16£ + 0.60 



R — J.-I-W.-T 




4 

s = = 250 

E,a n 




PO=lG0e-^ ¥ =2% 




P =-0.016- 0.012 j 




P l =-0.016 + 0.012 j 




p 2 =P 3 =P 4 =-0.16(or< 


-0.16) 



(27) 



Hence: 



The desired polynomial will be: 



(28) 



a(s) = (s + 0.016 + 0.012 j)(s + 0.016 - 0.012 j) 

(s + 0.2)(s + 0.3)(s + 0.4) 
We rewrite (29) as: 

a(s) = s 5 + 0.932s 4 + 0.2892s 3 + 0.0327s 2 + 

0.0009s + 0.0001 
Fro m (23), we write the vector gain K c as : 

K =[0.00001 0.0002 0.0121 0.1431 0.7407] 
Finally, the desired feedback vector is: 

K = KT 



(29) 



(30) 



(31) 



VI. Results of Simulation 

We perform some numerical simulations to illustrate the 
effectiveness and the robustness of the proposed methodology 
of control. 

In this section, we present some results using the given 
model for the naval vessel and the environmental disturbances. 

All the results are computed with the aid of MATLAB 
software. 

A. Directional stability 

If we look to the state matrix A of the considered system, 
the third row vector is null, this indicates that the system has a 
null eigenvalue. As result, the system is not stable. 

If the yaw angle is changed by an external force, the open- 
loop system changes the direction and continues moving on 
this direction; however in close-loop case the system returns to 
the initial direction. We can see the influence of the feedback 
on the trajectory in the Fig. 4. 

For the reader's convenience, the mathematical trajectory of 
the underactuated ship moving in surge, sway, and yaw (Fig. 4) 
is recaptured as: 



Xpos = (u cos(i//) - v sin(y/))dt 
Ypos = \(u sm(y/) - v cos(t//))dt 



(32) 



We can interpret the directional stability in the close-loop 
system by the asymptotic stability in yaw mode (Fig. 5). 
However the roll mode is stable in both cases, but the 
amplitude of damping increases in close-loop case (Fig. 6). 
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Figure 4 Directional stability of the trajectory 
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Figure 5 Yaw angle response in open and close-loop 
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Figure 6 roll angle response in open and closed loop 

Note that the initial conditions didn't present the 
equilibrium state (x(0) * 0). For this reason, the yaw in open- 
loop case changes the direction, however for the close-loop 
case the ship return to equilibrium state (see the trajectory in 
Fig. 4). 

For information, the initial state of the ship is taken as: 



v = -1.5, r = 0.01, y/- 



10 



p = 0.001, 



60 



with u„ 



B. Shaping yaw dynamics 

In addition to closed-loop asymptotic stability, we are often 
interested in other characteristics of the closed-loop transient 
response ( t , t p ,PO, and t). 

For instance, if we specify that the desired yaw response by 
PID controller and state-space feedback controller should have 
the smallest time constant with 2% overshoot, we get the 
results of Fig. 5. 

Note that the yaw reference is taken: y/ (t) = 100° . 

The initial state of the ship is: 

v = 0,r = 0,y/ = 0,p = 0,<j) = with u = 5 m.s 1 
It is clear in the yaw response (Fig. 7) that we can reduce 
the rise time from t R £l45s to t R = 110s by implementing a 
linear feedback controller instead of classical PID controller. 

The yaw response in the case PID controller presents also 
the criteria of minimal rise time t R with overshoot 
PO = 2% (Fig. 7 includes details). 

The parameters of the PID controller are: 

K p = -0.23, K, = 0.001, K D = 0.09 
However the amplitude of the roll damping increases for 
the feedback gain controller, which requires having a good 
equipment and efficient law to control roll motion (usually use 
fins to control the roll mode). 

Remember that the pair (A, B)is controllable, which is a 
necessary condition to apply the state-space feedback controller. 
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Figure 7 Shaping yaw response 
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Figure 8 Control influence on the rollmode 



C. Ship path following 

In this section, simulation results will be presented to prove 
the performances of the proposed feedback autopilot. For this 
purpose, we try to force the ship to follow a desired path 
(zigzag path) and try to apply wave disturbances on the motion. 

The results of simulation will prove the robustness of the 
two controllers in presence of disturbances. 

Briefly, the approach is to steer the ship such that it 
eliminates the distance between itself and the desired path (Fig. 
9). 

We define the following variables to mathematically 
formulate the control objective: 



y t =y t -y 



Where: 



W d 



■■ arcsin 



y„ 



H + yl 



(33) 



(34) 



The undesirable motion of a ship in a seaway is induced by 
the action of environmental disturbances: waves, wind and 
current. However, ocean waves are the dominant 
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environmental disturbances. As a result, the wave disturbances 
are considered at the output as it is shown in Fig. 10. 

Waves' disturbances can be characterized in terms of their 
frequency spectrum This frequency spectrum can be simulated 
by either a series of added sinusoidal signals, or using filtered 
white noise. 

In our work, the used filter to approximate the spectrum is 
of the form Re f. [11]: 



F(s) 



■0.1s + 0.25 



(35) 
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Figure 9 Deflection of the ship 
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Figure 10 System structure with disturbances 

Figs. 11 and 12 show X -Y plots of the ship zigzag 
trajectory; the test is done for the two different autopilots. 

For the test of tracking, we generate a zigzag path as 
desired trajectory. This test is a standard maneuver used to 
compare the maneuvering properties and control characteristics 
of a ship with those of other ships. Another feature is that the 
experimental results of the test can be used to calculate the 
parameters of Nomoto model. 
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Figure 11 Zigzag path following with and without low waves disturbance 
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Figure 12 Zigzag path following with high waves disturbance 

The results of two simulation tests are given here. The first 
simulation demonstrates the effectiveness of the feedback 
controller in absence of disturbance and its robustness 
comparatively with the PID autopilot in presence of 
disturbances (Fig. 11). In the second simulation, the influence 
of strong waves on the path following is demonstrated (Fig. 12). 

Fromthe first simulation results (Fig. 11) and in absence of 
disturbances, it can be concluded that feedback controller is 
more effective and faster than PID controller especially in 
terms of following the desired path. After adding small 
disturbances (white noise driving the filter has a power spectral 

density a 1 = 0.5), the feedback autopilot keep the following of 
the path, however the PID controller can only keep the path 
direction not the path (keep only the yaw direction). 

If we increase the waves' power ( a 2 =1) the influence on 
the ship was presented in Fig. 12. In this case, the robustness of 
the feedback gain controller is completely proved. 

VII. Conclusions 

The method for constructing ship autopilots based on linear 
feedback gain has been presented. The advantages of this 
method are to construct the controller easily and to have some 
malleability to shape the yaw dynamics as desired. Simulation 
results with real ship parameters show that by linear feedback 
controller we decrease the time constant comparatively with 
PID controller, which is an important characteristic especially 
for ship maneuvering. 

Simulation results show that this feedback autopilot is 
effective and robust in comparison with conventional autopilot 
based on PID control especially in the presence environmental 
disturbances (sea waves). 

The benefits of the proposed approach lie in the ability of 
such desired polynomial to better shape the yaw response as 
desired and to reject external disturbances. 

The presented approach can be readily applied to other 
vehicles or extended to higher dimensional control and 
guidance problems if the system is controllable. 

Among the problems of decreasing the yaw rise time, we 
mention the increasing of the roll damping amplitude; this 
problem is due to the roll coupling. And it can be settled by 
studying yaw response and roll controllers together. 
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Another limitation of the state-space feedback controller is 
the controllability of the system which is necessary. 
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